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Abstract 



A recursion technique of obtaining the asymptotical expansions for the 
bound-state energy eigenvalues of the radial Schrodinger equation with a 
position-dependent mass is presented. As an example of the application we 
\Q ■ calculate the energy eigenvalues for the Coulomb potential in the presence of 

\ position-dependent mass and we derive the inequalities regulating the shifts 

of the energy levels from their constant-mass positions. 
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X ■ 

Introduction 



The notion of effective mass, which was introduced more than 50 years ago to 
describe the motion of electrons in crystals pQ, has now a broader scope. It is 
employed in nuclear physics [2], theory of quantum liquids [3], metallic clusters [lj 
and semiconductor heterostructures: quantum wells, wires and dots [5j. Attempts 
to accommodate spatial inhomogeneity of multi-layer nanostructures have invoked 
treatment of the effective mass as a position-dependent quantity [fH E]. 

In order to calculate the bound-state spectrum of the Schrodinger equation with 
a position-dependent mass (PDM), various iterative [HI \10\, variational (TTJ [12] and 
perturbative [13] schemes were proposed. Usually, these schemes provide solutions in 
numerical form. On the other hand, there exist methods of deriving exact analytical 
solutions to the PDM Schrodinger equation for some potentials such as the Coulomb, 
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oscillator and Morse potentials [HI [TBJ [TB]. Having a solution in analytical form 
greatly simplifies investigation of the structure of energy levels. Thus an important 
task is the development of approximate analytical methods to solve the bound-state 
problem for the PDM Schrodinger equation. 

For the constant-mass Schrodinger equation with spherical symmetry, the 1/N- 
expansion method has proved to be efficient. The name of 1/iV-expansion refers to 
the group of related approaches [UJ HHJ [191 EH EI] that start with describing clas- 
sical motion of a particle located at the bottom of the effective-potential well and 
then consider quantum corrections. The versions of this method differ in choosing a 
formal expansion parameter. Namely, going to the classical limit, one takes for the 
small parameter the inverse number of spatial dimensions, the inverse principal quan- 
tum number, or the Planck constant. The last choice, the so-called ^-expansion [20J, 
seems to be the most natural one and, besides, allows to formalize calculations in a 
simple recursion procedure. 

The purpose of this work is to extend the /i-expansion technique to the PDM 
Schrodinger equation. We develop the procedure for calculating the bound-state 
energy eigenvalues and then apply it to investigate the influence of PDM on the 
Coulomb potential spectrum. 
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1 /i-expansion technique 

The PDM Hamiltonian of the particle is given by [22] 

h 2 

H = T + V(r), T = -— [m(r) a Vm(r) /3 Vm(r) 7 + m(r) 7 Vm(r) /3 Vm(r) a ] , (1) 

where m(r) is the PDM, V(r) is the interaction potential, a, /3, 7 are the ambiguity 
parameters (a + /3 + 7 = — 1). 

Let us consider the spherically symmetric case in which the stationary Schrodinger 
equation with the Hamiltonian ([1]) is reduced to the radial equation 

rf 2 + m'(r)/l d\ 1(1 + 1) 2m ( r )nj( r ) E)X?lj(r)-0 
dr 2 m(r) \r dr ) r 2 h 2 J 

i"(r) 2m'(r)\ , N m /2 (r)l . . 

m*\r) m z {r)r j m 6 [r) 

with I being the orbital quantum number and primes denoting derivatives with 
respect to r = |r|. 

Assuming V(r) and m(r) are analytical functions such that Eq. (J2]) possesses 
bound states, we intend to find the discrete energy spectrum. To that end, we 
generalize the ^-expansion technique [20] to the PDM case. 

First we recast Eq. ([2]) into the Riccati equation for the logarithmic derivative 
of the wave function. Upon substituting ip(r) = x{ r )\/ m ( r ) to eliminate the term 
with dip{r)/dr and then putting C{r) = fox'( r )/x( r ), we g e t 

hC\r) + C 2 (r) = + l) + (- A + a 7 + a + 7) h 2 Q 2 (r) - ^Q(r) - 



1 + a + 7 2 



H P(r) + 2m(r)(V(r) — E) (3) 
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where we denoted Q(r) = m'(r)/m(r), P(r) = (m"(r) + 2m'(r)/r)/m(r). 

We try to obtain the solution to the Riccati equation with series expansions in 
the Planck constant h 

oo oo 

C(r) = Ck(r)h k , E = J2 E k h k . (4) 

In view of the freedom in defining the classical limit, which is formally the limit 
of h — y 0, the centrifugal contribution can be rewritten as follows 

h 2 l{l + 1) = A 2 + hAA + h 2 B (5) 

that enables to consider different modifications of the method depending on values 
of parameters A and B [20]. A particular choice of these values amounts to fixing 
a zeroth approximation. In practice, one first specifies V(r) and m(r) and then 
choose A and B in such a way that higher-order corrections in the ^-expansions 
for eigenenergies (jlj) become as small as possible. Note that one may achieve yet 
faster convergence of this series expansion by means of the renormalization: rewrit- 
ing V(r) or m(r) as the fc-expansion and properly adjusting its coefficients. Such 
generalizations of the method lie however out of the scope of the present article. 

Substituting the expansions (j3J into the Riccati equation ([3]) and equating to 
zero coefficients of successive powers of h, we obtain the system of equations 

C 2 = 2m(r)(V(r)-E ) + ^, 



2 

i 



C + 2C\A = -2m{r)E x + 7l (^)' 
C[ + 2C C 2 + C\ = -2m(r)E 2 + j 2 (-) 2 + F(i 



' \ r 



2 

Cfc-i + Yl C * C k-* = -2™(r)E k + 7 fc (-) , k > 2. (6) 

i=0 



Here 71 = AA/r 2 , 72 = B/r 2 , 73 = 74 = . . . = 0, 



3 i ^ N \ s 1 + a + 7 



F(r) = I- + ai + a + 7 J Q\r) - -Q(r) '-P(r). (7) 

In order to take into account the nodes of the wave function for radially ex- 
cited states, we employ the Zwaan-Dunham quantization condition [231 124"] which 
expresses the principle of argument in complex analysis. For the discrete spectrum, 
solutions to Eq. (J2J) are real on the real axis and possess a finite number of simple 
nodes on it. The number of these nodes is just the radial quantum number n r . Then 
the integration of the logarithmic derivative C(r) along the contour that encloses 
only the above nodes yields 



1 

2tH 



C(r)dr = n r h, n r = 0, 1, 2, 



3 



This condition must be supplemented with a rule of achieving the classical limit 
h — > for the radial and orbital quantum numbers which are the specific quantum 
notions. We choose this rule as follows 



h — > 0, n r = const, I — > oo, hn r — > 0, hi = const. (9) 

Physically, it means that in the classical limit the particle is located in the minimum 
of the effective potential V(r) + A 2 /(2m(r)r 2 ) and thus moves along the stable 
circular orbit whose radius r is a solution to the equation 

m(r o y o V'(r ) = A 2 (l + . (10) 



Meanwhile, the energy in the zeroth approximation is given by 

A 2 

2m(ro)r 2 

and the quantization rule transforms into 

1 

2rd 



Eo = V(r ) + ^- 7Z — 2 (11) 



C k (r)dr = n r 5 k<1 . (12) 



With the help of such rules we can develop an algebraic recursion procedure for 
calculation of the ^.-expansion corrections to energies avoiding bulky formulas that 
arise in higher orders of the standard Rayleigh-Schrodinger perturbation theory. 



2 Recursion procedure for calculating energies 

To solve Eqs. ([6]), we expand the functions V(r), m(r), F(r), P{r) and Q(r) into 
Taylor series in powers of the new variable, the deviation from the minimum of the 
effective potential, x = (r — r )/ro, 



oo oo 



V(r) = s ^V i x\ m(r) = S ^m i x\ F(r) = s ^F i x\ 

i=0 i=0 i=0 

oo oo 

P(r) = ^P i x i , Q{r) = Y,QrX l (13) 



i=0 i=0 



where Vi = r % V^{rQ)/i\, mi = r l m^\r )/i\, whereas F^, Pi and Qi are expressed 
through Vi, rrii. 

Consider the first equation of Eqs. ([6]) that determines Co(r). Writing down the 
squared quantity as 



CjJ(r) = wV(l + clix + a 2 x 2 + . . .) (14) 



with 



u = \ l2(m V 2 + miVi) + (3m - m 2 ) 

2nio + mi 



Cli = 



2 - ( £ mj V,- H2 + ( (i + 3)(-l)' -^)), (15) 



we find the coefficients of the Taylor series expansion for Cq(x) 

oo 

Cq(x) = —ux\J 1 + a\X + a 2 x 2 + . . . = x C®x\ (16) 



i=0 

/- ! 

2^ 1 



CS = -", C? = ^(53^-0;^). (17) 



The " — " sign in Eq. (fT6l) is chosen so as to meet the boundary conditions: C (r) > 
for r < r , Cq(t) < for r > r , which follow from the quadratic integrability of the 
wave function. 

The function Cq(x) has a simple zero at x — 0. Then, after examining the system 
of equations (jEJ), we conclude that Ck(x) has a pole of the order of (2k — 1) at this 
point and, hence, can be represented by the Laurent series 



C k (x) =x 1 - 2k J2 C i xi - 



i=0 



With definition of residues, this expansion permits us to express the quantization 
conditions ( fT2l) explicitly in terms of the coefficients C\ 



C 2k-2 = —Sl,k- (19) 

Inserting the above expansions into Eqs. (|SJ) and collecting coefficients of powers 
of x, we obtain recursion formulas for calculating the yet undetermined C\ and also 
the corrections to eigenenergies Ek 



Cf 



2C° 



9(2 -2k + i)\ j k (~l)V ~2k + i) + E^ 2 6 K2 - 



k 2k-2 



m i+2 -2k 1 ^k-i S^S^^i^k-j 11 3-2k + i 



k—1 i i 
j=l p=0 p=l 



i y£ 2k - 2, (20) 



/ fc 2fc-2 \ 

^ = r- U + *b* W - -C£± 2 C i C 2^2- P (21) 

where B(/c) is the Heaviside step-like function (Q(k) = 1 for k > 0). 

The derived formulas completely resolve the task of calculating the energy spec- 
trum of the PDM Schrodinger equation. 

As an illustration, let us compute the energy eigenvalues in the field of the 
Coulomb potential V(r) = — q/r in the presence of the PDM 

m ( r ) = m ^ \ \ ( 22 ) 

v ' (1 + ar) x K ' 
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where the constants q, m c and a are assumed to be positive. We choose the paramet- 
ers A, A and B in the form 

A = h(n r + 1 + 1), A= -2n r -l, B = n 2 r + n r , (23) 

which in the constant-mass case guarantees that the Balmer formula for energies is 
reproduced just in the zeroth approximation Eq of the ^-expansion method. The 
ambiguity parameters are set to be a = 7 = 0, /3 = —1 in correspondence with the 
kinetic part of the effective Hamiltonian obtained for a crystal with slowly-varying 
inhomogeneity [6l [25] . 

The calculated partial sums E^ k ' = Eq + E\h+ . . . + Ekh k of the ^-expansions for 
energies as well as the exact energy values -E n um found by numerical integration are 
listed for various A in Table 1. The calculation has been carried out for the states 
with n r = 0, I = 2 and n r — 1, I — 1, using q — 10, a — 0.1 in units in which 
h = 2m c = 1. 



Table 1: Absolute values of partial sums E^ of /i-expansions for energies in the 
Coulomb potential in the presence of the PDM m(r) = m c /(l + ar) x . E num is the 
result of numerical integration. 
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,1 = 1) 




A = 2 
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A = -2 


A = -3 


A = 2 


A = 3 


A = -2 


A = -3 
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1.77778 


1.18817 


3.64395 


4.04566 


1.77778 


1.18817 


3.64395 


4.04566 


1 


1.94444 


1.45345 


3.50153 


3.83753 


2.27778 


1.98401 


3.21669 


3.42127 


2 


1.83333 


1.25876 


3.58047 


3.94991 


2.07556 


1.66974 


3.39891 


3.69143 


3 


1.83000 


1.24978 


3.57934 


3.94732 


2.05556 


1.61180 


3.39064 


3.67228 


4 


1.83111 


1.25190 


3.58015 


3.94898 


2.06000 


1.62647 


3.39352 


3.67887 


5 


1.83111 


1.25184 


3.58014 


3.94896 


2.06000 


1.62478 


3.39330 


3.67837 




1.83111 


1.25183 


3.58014 


3.94897 


mm 

2.06000 


1.62510 


3.39329 


3.67834 



It should be added that in the particular case of A = 2 the problem has an exact 
solution [13] and the ^-expansion technique restores it. Indeed, in this case we have 
£5 = Eq — . . . — and E^ coincides with the exact value, as seen from Table 1. 

3 Order of energy levels 

Let us apply the obtained recursion formulas to determine the order of energy levels 
in the Coulomb field in the presence of PDM. 

First of all, we notice that for the chosen values of the ambiguity parameters 
a = 7 = 0,/3 = — 1 there exists an inequality making comparison of the energy 
levels with their constant-mass positions. Namely, if m(r) and m c are the PDM and 
the constant-mass respectively, the corresponding energy eigenvalues E(n r , I) and 
E c (n r , I) satisfy 

E(n r ,l) ^ E c (n r ,l) if Vr 6 [0, 00) m(r) § m c , (24) 
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provided that these eigenvalues exist. 

The inequality (j24p holds valid not only for the Coulomb, but also for an arbitrary 
attractive potential since it is a direct consequence of the variational definition of the 
discrete spectrum of a Hamiltonian bounded below [26]. In particular, this inequality 
explains the increase of binding energy for impurity bound states in quantum dots 
that was predicted in [27] . 

To obtain further inequalities, we fix the shape of the PDM ( 1221) and consider 
the first and second ^.-expansion corrections whose contribution to the energy is 
decisive. According to Eqs. ffTTl) and fT2T|) . we have 



r 2m c rg 



where 



= (2n r + l)(l + ar ) x (u J r -A) 
1 2m c r 2 l 



m c q(l + ar )-^![a 2 r 2 (2 - A)(l - A) + 2ar (2 - A) + 2] 

r (2 + (2 - A)ar ) " (26) 



It is convenient to rewrite the sum of these two corrections in the form 
m c q 



E + hE 1 = E B + 



- r ,A 1 , / ^2 , ft(2n r + l)(u;ro/A-l) > 



2s 2 A 2 



t A - 1 + (s - If + v r ' ^ ' * 

A 



(27) 



with E-q = —m c q 2 /2h 2 {n r + I + l) 2 being the constant-mass energy, t — 1 + ar , 
s = t A -|(t A -t A - 1 ). 

Using the expression (1271) . we can establish the relative order of energy levels 
that refer to the same value of the principal quantum number n — n r + 1 + 1 and are 
degenerate in the constant-mass case. The quantities A, r , to, t and s, appeared in 
Eq. (|27|) . depend on n, but not on n r and I separately, so they cannot remove the 
degeneracy. It is however removed by the last term in the braces due to the factor 
(2n r + 1). The sign of this term coincides with that of the quantity (euro/ A — 1) for 
which we find 

^r _ u 2 r 2 /A 2 - 1 _ Xar [2 + (3 - A)ar ] 
A a;r /A + l 2(wr /A + l)(l + ar ) 2 ' 1 ' 

Obviously, if A < 0, then (ur /A — 1) > and, as a consequence, the level with 
the quantum numbers (n r , I) has the higher energy than the level with (n r — 1,1 + 1). 

Now let us prove that for A > the relative order of those levels is converse, i.e., 
E(n r ,l) < E(n r — 1,1 + 1). As seen from Eq. ( 1281) . it is sufficient to demonstrate 
that one has 2 + (3 — A)ar > when A > 0. Rewriting the equation ( |T0|) for r as 



A 



2 r 



m c q 



2 V 



A 



2 



m c q 



-s, (29) 



we see that it has a positive root only if s > 0. This entails ar Q = t — 1 < 2/(A — 2) 
for A > and, hence, 2 + (3 — A)ar > ar Q > that completes the proof. 
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Thus the inclusion of the PDM m(r) = m c /(l + ar) x removes the "accidental" 
degeneracy in the Coulomb spectrum so that 

E(n r ,l) ^E(n r - 1,1 + 1) if A ^ 0. (30) 

Interestingly, in the case of A < the above order of the energy levels coincides with 
that observed in muonic atoms and also in the soft-core Coulomb potential [28J. 

To proceed, let us consider spacings between the levels with the same value of 
n = n r + I + 1 and evaluate the ratio 

E(n r , I) - E(n r + 1, 1 - 1) ' 1 ' 

We adopt the expression for the energy up to the third correction, i.e., E = 
E + E\K + E 2 h 2 . Here E 2 calculated according to Eq. f l2T]) reads 

E * = 77 + n r) -8A 2 (ar (A - 2) - 2) 2 - 8Awr 6(1 + a x )- 

3a 2 r 2 (A - 2)(1 + A + ai) + ar (12(l + oi) - A(3ai - 4))) - 

w 2 r 2 (a 2 r 2 (16A 2 + 15a 2 + 8A(3ai + 1) - 12a 2 - 8) + 2ar (15a 2 + 12Aa x - 

12a 2 - 8) + 15a 2 - 12a 2 - 8)] + f{n) = + n r )g + f(n) (32) 

where f(n) is a function that depends on n, but not on n r and I separately. Its 
explicit form is inessential for us since we consider the levels with the same n. 
Combining Eqs. ( 1271) and ( |32l) . we find the requisite ratio 

1 + (2n r - l)6i/6g 
1 + (2n r + l)h/b 2 

where 

6 1 = 2^4' 34 

lDm c o; z rQ 

~ m c q 2 t x h / ur \ t A ~ 2 ^ 1 

A 3 s 2 V A 7 16m c w 2 r^ ' 1 J 

Clearly, one has i? < 1 for 61/62 > 0, R > 1 for 61/62 < 0. We managed to deduce 
the sign of 61/62 analytically in the case of the slowly- varying PDM. Assuming a is 
small and performing expansions in its powers, we get 61/62 = H/ (h — 2 A) + 0(a). 
Then from Eq. ([23D it follows that 6 x /6 2 < 0. 

Thus we conjecture the inequality 

E(n r l,l + l)-E(n r ,l) 
E(n r ,l)-E(n r + 1,1-1) ' 1 j 

It correctness was checked by an exact calculation of the energy via numerical in- 
tegration for various, not exclusively small, values of a and A. This inequality says 
that the shift on the energy scale is greater for levels with larger orbital quantum 
number. 

It is worthy to add that a similar problem, the PDM Dirac equation with the 
Coulomb potential, was considered in [291 [30]. There it was shown that relativistic 
factors such as an additional spin-orbit coupling due to the PDM also affect the 
order of energy levels. 



S 



Conclusion 



In the work we have developed the technique for the approximate calculation of the 
energy eigenvalues of the radial PDM Schrodinger equation. Based on the expansions 
in powers of the Planck constant with the subsequent analysis of the system near the 
effective-potential minimum, this technique provides, in principle, the calculation of 
the energy up to an arbitrary order in the analytical or numerical form. 

As an example, the spectrum of the Coulomb potential in the presence of the 
PDM m(r) = m c / (1 + ar) x has been investigated and the inequalities regulating the 
order of the energy levels in this system have been established. In particular, the 
inequalities ( 15Uj) and f[5Bj) show how the "accidental" degeneracy of the Coulomb 
spectrum is removed due to the PDM. 

The expressions for energies obtained with the developed technique can also 
be applied to study the dependence of the spectra of spherically symmetric nano- 
structures on their radiuses that is an actual problem of modern nanophysics. 
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